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In an earlier paper presented at ICES in 2015, we investigated solar particle event (SPE) 
radiation exposures (absorbed dose) to small, thinly-shielded spacecraft during a period when the 
monthly smoothed sunspot number (SSN) was less than 30. Although such months are generally 


considered “solar-quiet”, SPEs observed during these months even include Ground Level Events, 
the most energetic type of SPE. In this paper, we add to previous study those SPEs that occurred in 
1973-2015 when the SSN was greater than 30 but less than 50. Based on the observable energy range 
of the solar protons, we classify the event as GLEs, sub-GLEs, and sub-sub-GLEs, all of which are 
potential contributors to the radiation hazard. We use the spectra of these events to construct a 
probabilistic model of the absorbed dose due to solar protons when SSN < 50 at various confidence 
levels for various depths of shielding and for various mission durations. We provide plots and tables 
of solar proton-induced absorbed dose as functions of confidence level, shielding thickness, and 
mission-duration that will be useful to system designers. 
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ground level event (or enhancement) 
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unit of energy Mega-electron volt 
root mean square 

single-event effect 

solar energetic particle 

solar particle event 

smoothed sunspot number 

total ionizing dose 
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Introduction 


Aisne NASA missions have been with large spacecraft undergoing long durations. However, under 


recent budget restrictions, there has been a paradigm shift in the types of missions NASA has undertaken to those 
with smaller vehicles and shorter durations. Since many of these missions are also under constrained budgets, new, 
low-cost, and efficient methods are being employed on these missions to make them successful. One area to 
consider for cost-efficiency is the design of the vehicle for radiation protection. 


Furthermore, the recent trends in solar activity show a decrease in this activity from what has historically been 
normal in the manned spaceflight timeframe. The typical method of radiation analysis has been to consider the 
worst-case scenario for the radiation environment and then to design the vehicle to withstand the dose from that 
worst-case scenario. Thus, with the low-cost, short duration missions that are becoming more standard at NASA, 
using a historically worst-case method of analysis for total 1onizing dose may lead to an overdesigned vehicle that is 
more costly than necessary. 


In a previous paper’, we began the process of developing a new method for radiation analysis by investigating 
the solar-proton environment during months in the historical record when the monthly smoothed sunspot number 
(hereafter “SSN”’) was less than 30. In that study, we found that the worst-case solar particle event might not be the 
worst-case scenario for relatively thin spacecraft. In this paper, we extend the work to consider events in the 
historical record with sunspot numbers between 30 and 50. Additionally, we present a new probabilistic model of 
the absorbed dose due to solar protons at various confidence levels, depths of shielding, and mission durations for 
periods when SSN < 50. 


Background and Motivation 


The two most prevalent sources of primary radiation in deep space are galactic cosmic radiation (GCRs) and 
solar particle events (SPEs). Since SPEs tend to be a greater source of total ionizing dose for short duration 
missions, the focus of this work is on SPEs. Solar activity tends to follow an eleven-year cycle in which the sun has 
a higher frequency and more intense SPEs during solar maximum, and a reduced frequency and intensity of SPEs 
during solar minimum. These phases of solar activity are evident in the Figure 1, which shows the sunspot numbers 
for recent Solar Cycles and compares the relative rate of big monthly solar-proton fluences. Of the 568 months in 
the histogram, 38 have a very large accumulation of solar protons, with >10’ protons/cm* at >30 MeV. On the 
other hand, among the 242 months with SSN < 50, only four had such large proton fluences. Thus, the probability 
of a month with a very large solar proton fluence is roughly four times larger overall than it is during months with 
SSN <50O, 1.e., 6.7 + 1.0 % versus 1.7 + 0.8%. 


Methodology 


Thirty-nine solar particle events occurred during 1973-2015 with 30 < SSN <50. Of these events, six were 
Ground Level Events (GLEs), which are extremely energetic SPEs having proton energies above 400 MeV and 
producing secondary particles that can be observed with ground station neutron monitors on the surface of the Earth. 
Seven of the events are denoted as sub-GLE events, which produce measureable proton fluxes up to at least 300 
MeV but do not produce detectable levels of secondary atmospheric particles. The remaining 26 sub-sub GLEs are 
even less energetic with an observable increase in protons at energies greater than 30 MeV, but no observable proton 
flux above 300 MeV. Table 1 gives for each event the four Band fit parameters~° to the event-integrated integral 
spectrum in rigidity, J(>R). For detailed information on data sources and how we extracted the Band fits, refer to 
section III in our previous paper’. As footnoted in the table, J(>R) in a few of the smaller events is better fit as an 
exponential, a single power-law, or a so-called Ellison-Ramaty’ form, which is a single power-law times an 
exponential. For convenience, Table Al in the Appendix (shown at the end of the paper) provides the spectral fit 
parameters for events with SSN < 30. Table Al supercedes the list given in Reference 1, with improved fits for 
several events and additional sub-sub-GLEs not included in the original study. 
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We show a few plots of the events in the table above: two each for the GLEs (Figure 2), the sub-GLEs (Figure 
3), and the sub-sub-GLEs (Figure 4). The four band fit parameters are shown within each plot. In general, these fits 
are valid only for proton energies greater than 10 MeV, which is adequate since shielding makes lower-energy 
protons irrelevant to the radiation hazard for spacecraft electronics. With the large dynamic ranges in these plots, it 
is difficult to judge visually the goodness of fit. To quantify the goodness of fit, we note in the upper-right corner of 
each panel the root-mean square (rms) width of the each fit’s distribution of fractional residuals, that is, the values of 
(datum-fit)/datum. For all 90 events in the combined SSN < 50 study, the average rms width is 16.8%. Because of 
the inherently larger uncertainties in fluences extracted from neutron-monitors, the average rms width for the eight 
GLEs is 28.4%. 


Dose Calculations 


We made dose computations with HZETRN 2005, a high-energy particle transport/dose code developed at 
Langley Research Center®. We chose to use the 2005 version of HZETRN for this study due to its ease of use and 
faster runtimes. We used aluminum as the spacecraft material and considered thicknesses of 1.1, 2, 3, 4, 5, 10, 20, 
and 50 g/cm*. The minimum thickness value corresponds to the range of a 10 MeV proton in aluminum. The other 
values provide a range that will allow interpolations and weighted averages in assessing the performance of realistic 
shielding distributions. Since we are primarily concerned with non-human missions for this study, we chose a 
silicon detector to represent sensitive electronics that would be susceptible to total ionizing dose for these types of 
robotic missions. 


For each event, the input to HZETRN is the event-integrated differential fluence (in protons/em’-MeV), 
examples of which are shown in Figures 5 and 6 below. From the integral fluence in rigidity J(@>R), the differential 
fluence in kinetic energy is calculated as dJ/dE = (dJ/dR)/B, where R is rigidity, E is kinetic energy, and f is the 
proton speed in units of light-speed. 


The compiled results of absorbed dose (cGy [Si]) as a function of aluminum thickness are shown in Table 2. We 
also highlight the highest dose in green for each thickness represented. Doses for events with SSN < 30 are given in 
Table A2 in the Appendix. In general, we find the highest doses are produced by GLEs. Therefore, for the case 
where SSN < 50, the worst-case dose to a thinly shielded spacecraft would be a GLE-type event, based on the 
historical record. However, there is one sub-GLE (1984 Apr 25) that produces more dose under some shielding 
depths than three of the GLEs (1977 Sep 24, 1977 Sep 19, and 1973 Apr 29). As shown in Figure 5, the 1984 Apr 
25 sub-GLE has a higher fluence than the 1977 Sep 19 GLE until approximately 150 MeV and a higher fluence than 
the 1973 Apr 29 and 1977 Sep 24 GLEs until approximately 400 MeV". 


As shown in Figure 6, there are even some sub-sub-GLEs that produce higher doses under light shielding than 
some GLEs. The 1974 Sep 11 and 2005 Jan 16 sub-sub-GLEs have a higher differential fluence over the 1973 Apr 
29 and 1977 Sep 24 GLEs until approximately 60 MeV and 80 MeV, respectively. These higher fluences correlate 
to the dose results shown in Table 2. Thus, the smaller events cannot be completely discounted and should be 
included as part of the overall assessment. 


When we compare these results with our previous paper, we see that there are more GLEs in the months with 30 
< SSN < 50 (six GLEs in 70 months) than in months with SSN < 30 (two GLEs in 173 months). Moreover, the 
doses from GLEs with 30 < SSN < 50 are generally larger than those with SSN < 30. Higher sunspot number 
generally correlates with greater solar activity and intensity of events, and our results reflect this tendency. 


As noted in the introduction, the standard method of TID estimation has been to use the historical “worst-case 
SPE” to determine the potential worst-case dose exposure to the vehicle. This worst-case SPE is typically the 
October 1989 series of events, which contained three GLEs and one energetic storm particle (ESP) event when the 


' However, it may be worthwhile to note the following: except for the 1984 Apr 25 event, the solar-proton data for 
events from 1973 November to 1985 December were taken from instruments aboard the IMP-8 spacecraft. IMP-8 
had very large gaps in its data recovery during the 1984 Apr 25 event. The spectrum shown here is therefore based 
on “cleaned” data from the GOES-5 spacecraft, obtained from SPENVIS. Unlike data from subsequent GOES 
spacecraft, the GOES-5 data are not fully corrected for secondary-response contamination’. 
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CME-driven shock reached and washed-over Earth. Table 3 compares the doses for October 1989 to the worst case 
series of events for SSN < 50, January 2005. Using October 1989 potentially overestimates the dose exposure in a 
mission undertaken during relative solar quiet conditions, even for shielding thicker than 5.0 g/cm’, where the 
October 1989 dose is not dominated by the ESP event. This usage would in turn lead to overdesign of the shielding 
for the vehicle, increasing the mass and potential cost to the mission. 


Probabilistic Modeling of the Solar Proton Radiation Hazard During Periods of Low Solar 
Activity 


We now use these proton fluences and dose calculations to build a probabilistic space radiation hazard model. 
This model answers questions, such as “Given a mission of N months duration during low solar-activity, what is the 
probability of accumulating solar proton dose greater than a specified value?” As a technical point, individual events 
are not the appropriate starting point for a probabilistic model. This fact is easily illustrated by some simple 
considerations: in 243 months, we have identified 90 potentially significant proton events. This corresponds to a 
mean rate of w = 90/243 = 0.369 events/month. With this mean rate, Poisson statistics say that the probability of 
getting three or more events in one month is 0.62%. For a sample of 243 months, we should expect 1.55 such 
occurrences. However, our data show 11 months with three or more events. With wu = 1.55, the probability of 
observing 11 such months is (again, using Poisson statistics) 7.6 x 10°’. Multi-event months are more frequent than 
implied by these simple considerations because a given active region moving across the face of the Sun often 
produces more than one SEP event. 


Proper handling of these associations among events requires that we add together events that occur within a 
specified time interval. This interval should be long enough to group together all events associated with a particular 
active region, but not so long as to combine episodes likely to be independent. One month is an appropriate time 
scale for this purpose. As shown in Figure 7, there is no correlation between the solar proton production in a given 
month and the solar proton production in the preceding month, at least not during periods of low-solar activity. 


Thus, the first step in building a probabilistic model is to add together all proton fluences and doses associated 
with a given month. The 50 months with non-zero SEP quantities derived from our 90 events are listed in Table 4a 
and 4b for solar proton fluences and absorbed dose, respectively. 


For each quantity of interest, we next rank the monthly values in decreasing order. Our study is based on 243 
months with SSN < 50. For each month, we therefore evaluate the corresponding cumulative probability, namely 
(244-r)/244, where r = 1-50 is the rank. A month’s rank is not necessarily the same for all quantities of interest. The 
“month-size” and cumulative probabilities can then be analyzed as log-normal distributions.'”'' Error! Reference 
source not found. 8 shows the log-normal distribution for monthly solar proton fluence > 60 MeV (protons/cm’). 
The steepening in the distribution at low fluence values below 10° protons/cm” is likely due to inefficiency in 
pulling small events out of the GCR and instrumental background. The flattening at the largest fluences is likely due 
to the small sample of these very large events. In the middle range, the data give a reliable estimate of the log- 
normal distribution’s mean and standard deviations. These general features are true of all of our log-normal plots. 
Finally, since we start with spectral fits to each of the events in this study, we can generate log-normal fits for any 
energy threshold of interest. We have produced similar plots for >30 MeV, >100 MeV, >300 MeV, and >500 MeV. 


Previous studies of solar particle hazard probabilities have been performed in terms of particle fluence at various 
energies. However, since we have extracted the proton spectrum for each event, we can calculate the dose due to 
each event and make a direct probability analysis of dose under various depths of shielding. Figure 9 shows the log- 
normal distribution of the monthly accumulated solar-proton-dose under 3.0 g/cm’. The horizontal dashed line 
represents the dose under 3.0 g/cm* from GCRs (all species, not just protons) as calculated from the Badhwar- 
O’Neill GCR model’ for the 1977 solar-minimum GCR environment’. Figure 9 also shows that there is only a ~5% 
chance that the accumulated dose due to solar protons will exceed the dose due to GCRs. Similar plots were 
generated for aluminum shielding values of other thicknesses. The probabilities of GCR- and solar-proton- dose 
equivalence varied from 10% under 1.1 g/cm’ to less than 1% under 50 g/cm’. 


* It is estimated that the GCR dose in the most recent solar-minimum (~2009) was larger than 1977 by about 20%. 
(P. M. O’Neill, private communication, 2016) 
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Probability distributions like those shown in Figures 8 and 9 are sufficient to characterize missions lasting for a 
single month. But since 1973, there have been four eras with SSN < 50 lasting for 57, 44, 54, and 88 consecutive 
months, in 1973-1977, 1984-1987, 1993-1998, and 2004-2011, respectively. It is therefore potentially useful to 
extend this analysis to longer duration missions, particularly if the most recent extended period of SSN <50 
represents likely solar behavior for coming decades. Multi-month missions can be addressed by using random- 
sampling techniques analogous to those in previous studies by Feynman et al. 1993 and Tylka et al. 1997. 
Specifically, our procedure is as follows: 


e Fora mission lasting N months, make N independent random samplings from the log-normal distribution. 

e Add up these N values to obtain a mission-accumulated value. 

e Repeat the above two steps for a large number of trials and rank the results to get the corresponding 
cumulative probability distribution. In this work, 10’ trials were performed in order to guarantee a smooth 
distribution. 

e Repeat the above three steps for each quantity of interest. 


Figure 10 shows a sample of results from these simulations. The left panel shows the probability of solar-proton 
dose under 2.0 g/cm’ Al shielding for mission durations of 1, 2, 3, 5, and 7 years. The right panel shows the 
probability of solar-proton dose for a one-year mission under Al shielding depths ranging from 1.1 to 50.0 g/cm’. 


Finally, Table 5 shows the upper-limit on the solar-proton dose at confidence levels between 70% and 99% for 
various depths of shielding and mission durations. Also tabulated are the corresponding GCR doses based on the 
1977 solar minimum model. 


Comparison to the JPL91 Solar-Proton Fluence Model 


12 


Our methodology is inspired by that developed for solar proton fluences in the JPL91 Model”. In particular, 


e We employ log-normal distributions whose means and standard deviations are derived from 
observations; 
We use random-sampling techniques to simulate mission-accumulated quantities; and 

e We group correlated SEP events together. In our work, we use “months” for this grouping. JPL91 used 
episodes, many of which last for more than 20 days and hence are analogous to our “months” 


However, there are important differences in our approach: 


e For each of our events, we extracted the proton spectrum, allowing us to do probabilistic studies for 
arbitrary proton energy bins as well as doses. JPL91, on the other hand, was constrained to a handful of 
hardware-defined energy channels; 

e In developing the probability rankings, JPL91 counts only the events selected for the study. In our 
work, based on months, we also count the months in which no SEP significant events were found. 

e Asaresult, JPL91 includes an “events per year” parameter. In JPL91 simulations, the number of events 
in a trial is sampled from a Poisson distribution with this mean. In our simulations, on the other hand, 
there are always 12 monthly samplings in each year. 

e Most importantly, JPL91 included all SEP events, including those observed during periods of very high 
solar activity; in our work, we restrict our analysis to events observed during low-activity periods, with 
SSN < 50. 


The similarities between our procedures and those of JPL91 allow us to compare parameters of the data sets and 
the log-normal distributions. These comparisons are shown in Table 6 for two channels (>30 MeV and >60 MeV) 
that are featured in both models. We note these differences: 


e JPL91 has a higher event rate per month. In our case, this quantity is simply the number of months with 
significant SEPs (50) divided by the total number of months with SSN <50 (243). 
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The mean “event” size is larger by a factor of 10° in JPL91. 

e The largest “month” in JPL91 is October 1989, whose fluences are nearly an order of magnitude larger 
than those in the largest month in our survey, January 2005. 

e Most important, the slopes of the log-normal distributions (as given by the “Log;y Normal Sigma” 
values) are much larger in our analysis. 


This difference in slopes has important implications for the application of our model, as shown in Figure 11. 
These plots show the log-normal fits from JPL91 (red) and this work (blue) for proton fluences >30 MeV (left) and 
>60 MeV (right). For probabilities below ~99%, our fluence values fall below those of JPL91, as expected. For very 
rare fluence levels (>99%) our values are higher than those from JPL91. This is not reasonable. This problem 
presumably reflects a breakdown of the log-normal distribution at very large fluences, perhaps due to physical 
limitations on the SEP production. 


This problem also affects the mission-accumulated fluences, as illustrated for simulated missions of 1-year 
duration in Figure 12. For confidence levels above 99% (that is, probabilities below ~1%), our simulated fluence 
values exceed those from JPL91. This is not reasonable; but it is also not necessarily surprising, in that the historical 
record is insufficient to address such rare occurrences. Thus, we conclude that our model results for confidence 
levels above 99% CL are not reliable and should not be used. 


Discussion 


Sunspots, of course, do not cause solar particle events. We simply use them as a general indicator of the level of 
solar activity, which includes SPEs. As noted, the Sun appears to be moving into a less active regime, with fewer 
sunspots, fewer and smaller SPEs'*, and decreased coronal and interplanetary magnetic fields’. The latter is 
responsible for the higher GCR levels in the most recent solar-minimum. Shock-acceleration theory also suggests 
that the weaker fields play a role in the reduction of the SPE radiation-hazard'®. 


Our probabilistic calculations are based on the historic record of SPEs during times of low sunspot numbers, 
derived using data from 20 of the last 43 years of nearly continuous monitoring. One inherent assumption in using 
our results is that SPE production in future periods of low sunspot number will be no worse than what we have seen 
so far in the Space Age. This seems likely, but our current understanding of solar activity does not guarantee it. 


Utilizing these probabilistic assessments also requires faith in predictions of future solar activity. It may be safe 
to assume that current low-activity levels will continue for the next year or two. Beyond that, however, is 
speculative. We have included in our calculations mission durations of 3, 5, and 7 years simply because such 
calculations have been offered in previous SEP probabilistic models'”'*. Although the most recent era of SSN < 50 
persisted for 88 months (February 2004 - May 2012), whether these extended calculations will be useful depends on 
future advancements in forecasting solar activity. 


In this study, we have focused on a probabilistic model of solar-proton induced dose. But our method of 
extracting proton spectra offers other opportunities as well. For example, predictions of the total number of a certain 
type of proton-induced single-event-effects (SEEs) over a mission can also be addressed. To do so, it would be 
necessary to propagate our solar proton spectra through the shielding distribution, convolve the resulting spectra 
with the SEE production cross-section to get the number of SEEs (cf. Tylka et al. 1996), and then analyze those 
results with the same probabilistic procedures we have employed for dose. We note, however, that this work 
provides no guidance on peak SEE rates, which are also an important engineering consideration. 


Finally, as we also have discussed, the October 1989 event sequence has often been used as a plausible worst- 
case SPE environment. It may therefore be useful to note that the probabilistic dose predictions in this work suggest 
that there is a ~3-4% probability of encountering an October-1989-like dose level in a one-year mission at SSN < 50 
for shielding levels between 1.1 and 50 g/cm’. 
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Figure 1: Left: Monthly smoothed sunspot numbers since 1954 plotted versus the month of the Cycle (for Solar Cycles 
19-24. Right: Histogram of observed monthly solar-proton fluence above 30 MeV since 1967, with the red histogram 
representing all months and the gold representing months with SSN <50. 
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Figure 2: Band function fits to two GLEs. Note that data points at less than 0.13 GV (less than 10 MeV), which are 
generally irrelevant to dose calculations for spacecraft electronics, are not included in the fits. 
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Figure 3: Band function fits to two sub-GLEs. 
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Figure 4: Band function fits to sub-sub-GLEs. 
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Figure 5: Comparison of several GLEs to a sub-GLE that has greater differential fluence over certain parts of the 
spectrum. 
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Figure 6: Comparison of several GLEs to two sub-sub-GLEs that have greater differential fluence over certain parts 
of the spectrum. 
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SSN < 50, 1973-2015, 243 months 


cc = 0.058 
9g cc = -0.031 (excluding upper limits) 
® 


Monthly Solar Proton Fluence, >30 MeV, /cm42 
oO) 
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Previous Month Solar Proton Fluence, >30 MeV, /cm’2 


Figure 7: Correlation of monthly solar proton fluence at >30 MeV versus the same quantity in the preceding month. 
Months for which our search criteria found no significant SEP events are represented here as upper-limits at 1000 
protons/cm*. The correlation coefficients in the upper right indicate that there is no significant correlation, at least for 
months with SSN < 50. 
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Figure 8: Log-normal fit to the monthly accumulated solar-proton fluence at >60 MeV. 
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SSN < 50, 244 months, 1973-2015 
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Figure 9: Log-normal fit to the monthly accumulated dose under 3.0 g/cm? of aluminum shielding. 
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Figure 10: Simulations of accumulated solar-proton dose. Left: probability of dose under 2.0 g/cm” Al shielding for 
missions durations of 1, 2,3, 5 and 7 years. Right: Accumulated dose versus probability for a one-year mission under 
various depths of Al shielding between 1.1 and 50 g/cm”. 
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Figure 11: A comparison of the log-normal analysis for monthly solar proton fluences at >30 MeV (left) and >60 MeV 
(right) in this study (blue) and in the JPL91 model (red). 
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Figure 12: A comparison of probability vs. mission accumulated SEP fluence at >30 MeV (left) and >60 MeV (right) in a 
1-year mission from this work (blue) and from the JPL 91 model (red). 
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Table 1: Spectral-Fit Parameters for Solar Proton Events 


during Months with 30 < SSN < 50 





| ces | | | | | ft jf 
| 1997 | Nov | 06 |35.0/ 815E+08 | 0284 | 538 | 01159 
34.6 

was | | | | | | | 
31.8 
49.8 
35.0 
35.3 
apse | [ | [| | [| [| 
34.4 
34.0 
32.1 
| 1977 | Sep | 08 {39.1 
39.1 
47.5 
46.5 
46.7 
Apr 
40.2 
Sep 
Nov 

Alternate functional forms: a=exponential; b=Ellison-Ramaty form; | 
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Table 2: Calculated Absorbed Doses for Solar Proton Events 


| rl CU Absorbed Dose (cGy) / Thickness of Aluminum Shielding (g/cn’) 

| Year =(|Mon/Day 14 | 2 | 3 | 4 | 5 | 10 | 2 | 50 
| Ges | | | {| | ff 
| 1997 | Nov | 06 | 6.18E+01) 3.33E+01| 2.11E+01] 1.49€+01] 1.12E+01 | 4.26E+00| 1.28E+00) 1.856-01 | 
|subGles | | [| | | | | | {| Jj | 
sub-sub-GLEs| | | | | 
| 1977 | Sep | 8 | 5.07E-01 | 1.32E-01 | 5.27E-02 | 2.76E-02 | 1.67E-02 | 3.61E-03 | 7.19E-04 | 8.72E-05 
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Table 3: Accumulated Solar-Proton Dose in October 1989 and January 2005 


| 
Event Type 
1989 | 2140 | 1220 | 816 | 604 | 473 | 209 | 784 | 154 | 
1989 
1989 
1989 
1989 
1S 0.25 
| 0.096 


| 8622 @8| 14 #| 0.618 0.096 0.0136 0.0013 
73 0.957 0.102 


ff 
2005 0.094 0.047 0.029 0.020 0.0063 0.0019 0.0004 
195.0 


m 1 
9 : 


NIN IN IN 
OQ;/O/O;/O 
O}O/O!}O 
nya] on} an 


| 5.5 | A 
1140 | 508 | 66.2 | 338 | 269 | 127 | 5.21 | 124 
327.2 130.3 105.7 | 564 | 421 | 169 | 617 | 134 | 
Ratio. October 1989 / January 2005: 5.97 | 674 | 467 | 576 | 556 | 497 | 436 | 370 | 


Table 4a: Monthly-Accumulated Solar Proton Fluences 
During Months with SSN < 50 


Monthly Solar Proton Fluence, cm*-2 


1973| Sep 34.4 | 3.04E+06| 7.25E+04| 1.55E+03| 2.33E-03 | 8.07E-08 
3.65E+04 
1974| Jul 34.0 | 8.62E+06|3.77E+05| 3.65E+04 | 2.05E+02 | 1.59E+01 


1974 Nov ar S 4A6E+06 | 6-47E+05 | 1.95E*05 | S Sle tOs | 6.50E+02 
LerfS Aug 14.3 24 OG | ASE tOS | L.O2ZE+OS | 2 OSE+03:| S.12b+O02 














N) 





3aiN) 


N) 





1975| Nov | 16.1 |2.13E+04|3.03E+03| 7.04E+02 | 2.63E+01 | 5.06E+00 


1376 Mar 12 .63E+05 | 2.51E+05| 1.47E+05 | 4.88E+04 | 2.83E+04 
LeTs Apr V2 .ffE+O#/ | 9. 24E+06 |] 3.54E+06 | 2.18E+05 | 3.48E+04 


NO 


— 


1976| Aug | 14.0 | 1.03E+06]| 2.12E+05| 4.77E+04 | 4.19E+02 | 2.80E+01 





3] )|N 





—/ (| © 


1985| Jan | 205 |1.19E+06]| 2.79E+05| 7.08E+04 | 8.30E+02 | 5.22E+01 


1985 Apr 13 ee | 4700S | *_26FtOs | S.35E70s: | 672 +O2 
1985 Jul 17.4 .89E+06 | 1.36E+06 |] 6.25E+05 | 1.45E+05 | 8.42E+04 

















ie 


1986 47E+07 | 7.70E+06| 2.22E+06 | 6.87E+04 | 6.87E+03 
1986| Mar 13.0 O8E+05 | 4.68E+04]/ 8.73E+03 | 1.75E+02 | 2.47E+01 
1986| May 14.3 64E+05 | 2.86E+05]/ 1.03E+05 | 1.62E+04 | 8.02E+03 


7 7 
76E+04 
| 
1994| Feb 34.8 | | | : | 

E04 
1994| Sep 26.6 : | 2.56E+01 | 4.10E+00 
.38E+06 
82E+03 
-15E+06 
6 S0E+03 
1997| Sep 28.3 |4.90E+04/2.77E+04| 1.81E+04 | 6 90E+03 | 4. 26E+03 


2004 6.37E+06 | 4.27E+05| 7.29E+04 | 3.63E+03 | 1.00E+03 
2004 5.37E+06 | 4.09E+05|/ 7.21E+04 | 1.84E+03 | 2.94E+02 


2004 Nov Soe 4. 50E+07 6.04E+06 1. 18E+06 SOGE+04.| SO2ZETGs 
39 | 2004 Dec as 1.20E+05 | 1./5E+04] 4.04E+03 | 1.50E+02 | 2.88E+01 





SIOIN 
\)| 
NO 
mM) mM 
+| + 
O/O 
OO 
N 
) 
O 
M7 
+ 
O 
On 
On 
n 
On 
M7 
+ 
O 
K 
W) 
O 
K 
M7 
+ 
O 
Q) 
| 
com 
oaks 
M7 
+ 
O 
N 


7eE+05 
1977| Oct 45.6 oe BS TS See TOd 1.43E+04 | 5.70E+02 | 1.13E+02 
1984] Jun 46.5 24E+04 | 5.91E+03] 1.61E+03 | 8.60E+01 | 1.98E+01 


a 
sap 





AIN 
W 
os 
M7 
iz 
O 
4 
6) 
O 
O 
MT 
x 
© 
W 
(0 
0 
2) 
M7 
+ 
© 
N 


2 
= 
2 
2 
2 
2 
2 
2 
= 
2 


i 


—/% 


N) 











2005 May 28.9 4.32E+06 | 347E+05 | 3.87E+04 | ..37E+0s3:| 203E+02 
2005 Jul 29.1 2.26E+OF | 1.66E+065 | 1.826705 | 1.52E703 
4/f | 2006 Dec A 4 4S8E+06 | 1.3TE+08 | 4.61TEt07 | 2: 37E+06 | 3.7 1TE+05 
oO: | 2011 Mar 36.9 4 19E*06 |S S56E+05 | “2605 | S32EF0s-| 1.146702 
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Table 4b: Monthly-Accumulated Absorbed Dose due to Solar Protons 
during Months with SSN < 50 


| _|Year| Month | __—__—__—Dose from Solar Protons, cGy (Si), under depths of Aluminum shielding 
| | | Atgicm2 | 2gicm2 | 3 gicm2 | 4gicm2 | Sglcm2 | 10 gicm2 | 20 gicm2 | 50 gicm2_ 


| 6 [1974] Nov | 1.20E+00 | 5.17E-01 | 2845-01 | 1836-01 | 1.28E-01 | 4.00E-02 | 1.01E-02 | 1.16E-03 _ 
| 8 [1975] Nov | 1376-02 | 4.03E-03 | 175E-03 | 9.72E-04 | 6.16E-04 | 1.53E-04 | 3495-05 | 4.88E-06_ 
| 9 | 1976] Mar | 2.13E-01 | 1.05E-01 | 655E-02 | 4.77E-02 | 377E-02 | 191E-02 | 994E-03 | 4.27E-03 
2.52E-06 
1.00E-05 
2.56E-04 
1986 1.20E-03 
1.30E-02 
1986 7.08E+00 1.45E-02 
6. 17E-0 4.04E-05 
2.03E-0 1.69E-03 
87 4.96E-0 5.64E-07 


co 


— 
© 


— | — 
}@© 


— 
(© 


J 
7.91E- 


5.52E-06 
Oct 2.79E- 6.42E-04 
Apr 1.83E-0 1.73E-06 
1995 1.87E-04 
3.96E-06 
2006 3.93E-04 
1.61E-01 
1.33E-04 
006/ Dec 4.21E-01 
8.33E-05 
1.46E-04 
2.72E+00 9.35E-04 


1.19E- 1.92E-04 
5 STE: 5 87E-04 


) ; 
0 ’ 

4.07E-0 3.99E-07 
0 : 
0 : 


—s: 
© 
(© 
& 


2 
| 
3 
| 
| 
3 
3 
| 
3 


— 


O89) DO) RO] ROPERS] ROT RST RO] NOT ROT RO | aS] mS] 
O| ©] OO} SI] | 1] B] OO] NO] =| O |] OO] oo] ~ 
—s 


NO 


Or; S]_ |e) || 
NO 
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Table 5: Model Parameters for JPL91 and This Work 
530 MeV LS OMe 


a 
Largest Largest 
Average Rate, Mined Event in 0 at Event in 
Events/month} ‘\2/'" Sample, ermal Sample, 


Sj FE / tt Sj 


| 060 | 14.10 |1.00E+07/5.02E+09// 039 | 1.07 |8.00E+06|2.69E+09 
This Work{/} 021 | 230 | 1.84E+03/7.21E+08|| 021 | 238 |260E+02|234E+08 


Table 6: Upper Limits on Solar-Proton-Induced Dose (cGy) in Years with SSN < 50 
(Values in red exceed the GCR-Induced Dose) 


2.59E+02 1.23E+03 
2.46E+03 
1.89E+03 7.25E+03 
1.97E+04 


5.67E+04 1.08E+05 1.68E+05 
48.96 114.24 


=e 
9.94E+01 2.64E+02 4,.86E+02 


2.216402 | —5.59E+02 | 9. 94E +02 
7.55E+02 1.78E+03 3.02E+03 
2.31E+03 5.17E+03 8.47E+03 


Depth of Confidence 
Shielding (Al Level 

1.1 g/em2 70% 
80% 
90% 
95% 
99% 


2 Years 


1.16E+02 
2.65E+02 
9.40E+02 
2.95E+03 
3.23E+04 
32.64 


1 Year 


2.65E+01 
6.73E+01 
2.70E+02 


1.17E+04 
16.32 


2.0 g/cm2 70% 
80% 
90% 
95% 
99% 
GCR 


9.74E+00 
2.53E+01 
1.05E+02 
3.75E+02 
4.94E+03 
16.32 


4.40E+01 
1.03E+02 
3.76E+02 
1.21E+03 
1.40E+04 
32.64 


2.43E+04 7.68E+04 
48.96 114.24 


es 
3.0 g/cem2 70% 5.80E+01 3.18E+02 
80% 
90% 
95% 
99% 
GCR 


4.86E+00 
1.37E+01 
6.36E+01 
2.50E+02 


4.03E+03 
16.32 


2.42E+01 
6.08E+01 
2.49E+02 
8.76E+02 


1.20E+04 
32.64 


1.40E+02 3.74E+02 6.95E+02 


1.78E+03 4.16E+03 7.O9E+03 
2.24E+04 4.73E+04 7.58E+04 


5.30E+02 1.31E+03 2.33E+03 


4.0 g/cm2 70% 
80% 
90% 
95% 
99% 
GCR 


2.83E+00 
8.05E+00 
3.81E+01 
1.53E+02 
2.53E+03 
16.32 


1.43E+01 
3.64E+01 
1.51E+02 
5.41E+02 
7 .62E+03 
32.64 


3.46E+01 1.94E+02 
8.43E+01 | =o 2.28E+02— 4.27E+02 
3.25E+02 1.46E+03 
1.10E+03 4.48E+03 
1.44E+04 4.98E+04 


70% 
80% 
90% 
95% 
99% 
GCR 


1.91E+00 
5.53E+00 
2.69E+01 
1.10E+02 
1.92E+03 
16.32 


9.90E+00 
2.56E+01 
1.09E+02 
3.99E+02 
5.88E+03 
32.64 


pO 
2.45E+01 7.13E+01 1.39E+02 


6.03E+01 3.12E+02 
2.39E+02 1.08E+03 
8.2 5E+02 3.40E+03 
1.12E+04 3.88E+04 


5.0 g/cm2 


48.96 114.24 


48.96 114.24 


Eee 
7.99E+00 4.99E+01 
2.07E+01 1.17E+02 
8.90E+01 4,.39E+02 
1.48E+03 


10 g/cem2 70% 
80% 
90% 
95% 
99% 
GCR 


5.50E-01 
1.70E+00 
9.02E+00 
3.99E+01 
8.11E+02 
16.08 


3.11E+00 
8.56E+00 
3.92E+01 
1.56E+02 
2.66E+03 
32.16 


5.12E+03 1.16E+04 1.98E+04 
48.24 112.56 


48.96 114.24 


eee 
2.83E+00 9.11E+00 1.92E+01 


7.73E+00 4.74E+01 
3.59E+01 1.92E+02 
1.42E+02 6.85E+02 
2.56E+03 1.01E+04 
46.44 108.36 
ee 
| 1.75E+00 


1.75E+00 


20 g/cm2 70% 
80% 
90% 
95% 
99% 
GCR 


1.72E-01 
5.62E-01 
3.26E+00 
1.56E+01 
3.68E+02 
15.48 


1.05E+00 
3.06E+00 
1.52E+01 
6.48E+01 
1.28E+03 
30.96 
50 g/cm2 70% 
80% 
90% 
95% 
99% 
GCR 


2.82E-02 
9.71E-02 
6.08E-01 
3.11E+00 
8.40E+01 
14.16 


1.84E-01 
5.63E-01 
3.00E+00 
1.37E+01 
3.01E+02 
28.32 
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5.14E-01 


3.82E+00 
1.48E+00 9.87E+00 
7.36E+00 4.23E+01 
3.12E+01 1.61E+02 
6.35E+02 2.68E+03 


42.48 99.12 





Ee —— 

co pe leone eee neem seen edema a | 
po CUMGCRU| 628A BO | 114.24 
I ae ed 
po CCGCRT | 16.2 8A 4B | 114.24 
ee | ee ee Oe a | 
| Cd UUGORLUTC G2 | 4 BS , 
Leeper seer eens eens seenenenepeeseereseerenneen | seennesersseeereerseeen a | 
po CUMGCORT UT GBA] 8A BOG | 114.24 
i Dd 
| UR ||| | ESS 
a ee | 

|  GCR | 15.48 | 46.44 | 108.36 

rr eee eatery 
| 60 giem2_ | 70% | 2 B2E-02 || 1.84E-01 | 5. 14E-01 | 3.82E+00_ 
po 8H | 97IE-O2 | 5. GBE-O1 | 1.48E+00 | 9.87E+00 | 
po 90% | 6 OBE-01 | 3.00E+00 | 7. 36E+00 | 4.235401 
po 9% SIEHOO | 1.37E+01 | 3.12E+01 | 1.G1E+02 
Pp 99% | BB AOEFO1 | S.O1E+O2 | 6. 35E+02 | os 2.68E+O3 
P| GOR 4S 28.32 424 Le 


Appendix 


Table Al contains spectral fit parameters for solar proton events that occurred when SSN < 30. This is an 
update of the results from Reference 1, with some revised fit parameters and additional sub-sub-GLEs that were not 
included in the original listing. Table A2 contains the corresponding calculated doses for these events. 


Table A1: Spectral-Fit Parameters for Solar Proton Events 


ng Months with SSN < 30 


dur 





























a a a ee ee ee ee 
| osub-Gtes {|| 
Aug 
| 1985 | ul || 09 [17.4] 236E+06 | 1.506 | 7.80 | 0.1140 
| 1986 | Feb | 06 [13.1]  1.16E+0O7 | 1.030 | 10.20 | 0.1425 
1986 Feb | 14 | 13.1 4.70E+05 
Sere | cag | oa as | Sees —| oer | waa} es — 
1975 Aug 04 |14.3 2.99E+06 4.44 
Nov 
n/a 
| 1986 | Feb | 04 /13.1] 1.29F+06 | 0684 | 533 | 0.0847 
7.37 0.0697 
| 1986] Mar | 06 [13.0] 1.41E+05 | 2.296 | 654 | 00949 © 
| 1985 | May | 08 | 14.3 n/a n/a 
°14987 May 29 | 26.5 1.62E+06 1.20 n/a 0.0427 
#14994 Jul 12 |28.5 2.67E+07 0.000 n/a 0.0340 
Oct 
0.447 7.67 0.0443 
1.980 n/a 0.1160 
n/a n/a 
2005 Jul 13 |29.1 2.85E+07 0.790 7.13 0.0404 
| C2005] ut || 14 [29.1] 5 96E+09 | 0.013 | 969 | 00374 
Sep n/a 
2006 4.92 0.0771 
2006 Dec | 14 | 12.1 4.67 0.0674 





Alternate functional forms:a 


exponential; b=Ellison-Ramaty form; c=>powrer-law 
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Table A2: Calculated Absorbed Doses for Solar Proton Events 


during Months with SSN < 30 


ff | Absorbed Dose (cGy) / Thickness of Aluminum Shielding (g/cm*) 

| Year |Mon/Day) 4.4 | 2) | 3 | 4 | 
ame | | | $=jW~ jf 7; ij ij Jf #7 i£=JJ{ | 

| subGlLEs | {| | | 

| 1985 | Jul | 9 | 1.27E+00| 5.34E-01 | 2.85E-01 | 1.79E-01 | 1.23E-01 | 3.54E-02 | 7.93E-03 | 7.39E-04| 

| 1986 | Feb| 6 | 4.88E+01] 2.43E+01| 1.47E+01| 1.02E+01 | 7.62E+00| 2.91E+00| 9.42E-01 | 1.24E-02 | 





| 2006 | Dec| 6 | 1.46E+02| 5.94E+01 | 3.15E+01| 1.98E+01 | 1.37E+01 | 4.18E+00] 1.06E+00| 1.28E-01 | 
jsub-sub-GLEs| | [| | |] ff 
| 1986 | Mar| 6 | 1.83E-01 | 6.17E-02| 2.80E-02 | 1.57E-02 | 9.76E-03 | 2.06E-03 | 3.74E-04 | 4.04E-05. 
| 1986 | May| 8 | 3.82E-02| 2.14E-02| 1.45E-02 | 1.11E-02 | 8.98E-03 | 4.79E-03 | 2.54E-03 | 1.09E-03. 
8.37E+00 
1.17E+00 
6.04E+00 
1.26E+01 
1.43E+00 
| 2006 | _ Jul | 6 | 1.55E-01 | 6.07E-02| 2.99E-02 | 1.75E-02 | 1.15E-02 | 3.20E-03 | 8.23E-04 | 1.33E-04 
‘Aug | 14 : 3 ; ; ; ; 


1.84E-01 | 5.98E-02 | 2.65E-02 | 1.49E-02 | 9.56E-03 | 2.45E-03 | 5./6E-04 | 8.33E-05 
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